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Abstract 

A Fermi resonance-algebraic model is proposed for molecular vibrations, where 
a U(2) algebra is used for describing the vibrations of each bond, and Fermi reso- 
nances between stretching and bending modes are taken into account. The model 
for a bent molecule XY2 and a molecule XY3 is successfully applied to fit the 
recently observed vibrational spectrum of the water molecule and arsine (AsHa), 
respectively, and results are compared with those of other models. Calculations 
show that algebraic approaches can be used as an effective method for describing 
molecular vibrations with small standard deviations. 
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Algebraic methods have been apphed to molecular systems for a number of years, 
and a recent review is available. U(4) and U(2) algebraic model have mostly been 
used so far for a description of rovibrational spectra of molecules. U(4) model has two 
advantages that rotations and vibrations are treated simultaneously, and that Fermi 
resonances are described by the nondiagonal matrix elements of Major ana operators. 
But it becomes quite complicated when the number of atoms in a molecule is larger 
than 4. U(2) model is particularly well suited for dealing with the stretching vibrations 
of polyatomic molecules, and currently under further development 

Recently, an important development on these methods is to incorporate bending 
vibrational modes of polyatomic molecules. Sample molecules such as acetylene '^1 , as 
well as X3 and XY4 [^'^1 molecules have been studied via algebraic methods. These 
approaches suffer from the neglect of the Fermi resonance interactions between stretch- 
ing and bending modes. Due to their importance as a mechanism for intramolecular 
energy transfer and for descriptions of excited stretching and bending vibrations, Fermi 
resonances have been taken into account in simple Fermi resonance-local mode models 
for bent triatomic and pyramidal XH3 1^1 systems and in boson-realization models 
for bent triatomic ^^'^1 and XY4 '^^1 molecules. In this paper we will introduce a Fermi 
resonance- algebraic model, where a U(2) algebra is used for describing the vibration 
of each bond. Our method is addressed for a bent molecule XY2 and a molecule XY3, 
where Fermi resonances between the stretch and the bend are considered. In a limit, 
the model corresponds to the boson-realization model. Vibrational spectra of two sam- 
ple molecules H2O and ASH3 are used to test our model. This method can be extended 
for describing the vibrations of molecules with other symmetry. The corresponding 
results will be discussed elsewhere. 



1. Hamiltonian 

Let n be the total degrees of freedom of vibrations and m the degrees of freedom 
of stretching vibrations in a polyatomic molecule. We introduce n U(2) algebras to 
describe vibrations: Uj(2) (1 < i < m) for stretching modes and U^(2) {1 + m < 
fjL < n) for bending modes. Each Ua(2) (1 < a < n) is generated by the operators 
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{Na, j+,aj-,a, jo,a}, Satisfying the following commutation relations: 




(1.1) 



where Na is related with the Casimir operator of U(2): 



2JI^ + J+,aJ-,a + J-,J+,a = M^a/^ + 1). 



The local basis states for each bond are written as \Na, Va), where Na is the eigen- 



Thosc of equivalent bonds are equal to each other: Ni = Ns, and N^ = N^, where 
footnotes s and b refer to the stretching and the bending vibrations, respectively. 

Due to Fermi resonances in the considered system, we will take V — Vs + Vb /2 as 
a conservative quantity, where Vs and Vf, is the total number of phonons of stretching 
and bending modes, respectively. Define the scale transformations for convenience: 



In the limit of A^s — > oo and A^b — > oo, 00(^0) reduce to bosonic operators with usual 
boson commutation relations PI. 

In the following, we will present the vibrational model in terms of Eq. (1.3) for two 
molecular systems with different symmetries. 

1.1 The model for vibrations of a bent molecule XY2 

For a bent triatomic molecule, we have n = 3 and m — 2. The Hamiltonian, where 
the interactions are taken up to the fifth orders, is given as follows: 



value of Na and Va the number of quanta on the ath bond. The action of J±^a on the 
local states is given by 




(1.2) 




(1.3) 
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H ^ alai{ujs + Xgaiol) + alasiujf, + Xbttaal) + Xi{a\a2 + 4ai) 

i=l 

+ \2{aiaia2a2) + \3{aiala2a2 + H.c.) + \4{a\a2 + a2a\){a\ai + 0202) 
+ X^{a\a2 + 0201)0303 + A6(o|oi + 0202)0303 + A7{(o| + 02)0303 + H.c.} 
+ A8{o|o2(oi + 02)0303 + if.c.} + A9{(o|o|oi + 020202)0303 + if.c.} 
+ Aio{(alala2 + 020301)0303 + H.c.} + AnKoj + 02)03030303 + H.c.}, 

(1.4) 

where uj and x are used for comparison with the constants in the Morse potential. 
Ai (1 < i < 11) arc the couphng constants. The term with A7 describes the Fermi 
resonance between the stretching and the bending vibrations. The last four terms in 
Eq. (1.4) are the fifth order interactions caused by Fermi resonance. In the limit of 
— > 00 and Nft — > 00, Eq. (1.4) corresponds to the boson-realization model '^^l 

1.2 The model for vibrations of a molecule XY3 

For a C^y symmetric molecule, we have n = 6 and m = 3. From the standard 
method of group theory, it is not difficult to obtain the character table, the represen- 
tation matrices of the generator, and the Clebsch-Gordan coefficients of the Ca^, group. 
From that knowledge, Hamiltonian where the interactions are taken up to the fourth 
orders are expressed in terms of operators Eq. (1.3): 

3 6 
H ^ a\ai{ujs + Xs Ojo|) + Y + ^6 S^iL) + ^ 4% 

i=l ^l=A i^j 

+ ^2 XI "m*^^ + ^3 (0I04O4 + a^a^a^ + oJoeOe + H.c.) 

+ ^74 {0.1(0505 -|- oeOe) + 02(0404 -|- ogOg) + 03(0404 -|- 05O5) -|- H.c.} 

7^5 {a|a4(a5 + Og) + 0305(04 + Og) 0303(04 + 05) H.c.} 
-I- 7^6 (01O5O6 + O2O4O6 + O3O4O5 -I- i^.c.) + ?77 X a}^ai(i'ja^ 

+ ^8 (o-^Oj -I- (^^a^^{a'^a^ -\- o]oj) + % (^^{(^{{o^^a}^ + o|.Oj) 

+ ?7io Y i4(4<^j<^j + H.c.) ?7ii Y + H.c.) 

(1.5) 

where lu and x are again used for comparison with the constants in the Morse potential. 
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7]i and 7/2 is the harmonic constant for the stretch and the bend, respectively. These r]p 
{3 < p < 6) arc the Fermi resonance couphng constants, Tjg {7 < q < 11) anharmonic 
constants of the stretching vibrations, and r}r (12 < r < 15) anharmonic constants 
of interactions between the stretch and the bend. In Eq. (1.5), we have omitted the 
interaction terms for the bending vibrations, which are similar to the terms with rjg 
(7 < g < 11) for the stretching modes, and also omitted the fourth-order interactions 
between the stretching and the bending vibrations. In the limit of — >■ oo and N;, 
— > oo, equation (1.5) reduces to the extended local mode model for three equivalent 
stretching bonds '^^1 , where the bending vibrations were neglected. 

2. Applications 

2.1 Application to H2O 

Studies of the vibrations of the water molecule already exist in the literature. Its 
stretching vibrational energy levels were calculated via nonlinear quantum theory l^^l . 
Halonen and Carrington proposed a simple vibrational curvilinear internal coordinate 
hamiltonian '^1 for its vibrational levels up to 18,500 cm^^. The same levels were 
studied by U(4) algebraic model '^l Its recently observed vibrational spectrum up to 
22,000 cm~^ were analyzed by the potential-energy surface from ab initio calculation 
l^^l as well as in the boson-realization model and the corresponding g-deformed model 
Owing to much research on this molecule, it provides a good testing ground for 
different models. 

Now, we calculate the Hamiltonian matrix elements in the symmetrized bases, then 
determine the 15 parameters in Eq. (1.4) by a least-squares optimization in fitting 
the experimental data Eo^g in cm~^ (their normal labels (i^i^'2i^3)), taken from the 
compilation of Ref. [14]. Two boson numbers and are taken to be 48 and 65, 
respectively. The parameters obtained are given in cm~^ as follows: 

ujs = 3705.121, Xs = -4.602, ujb = 1599.485, xi, = 3.913, Ai = -51.335, 
A2 = -12.746, A3 = -1.004, A4 = 3.199, A5 = -2.853, Ag = -21.956, 
A7 = 12.686, As = -0.422, A9 = -5.901, Aio = 2.187, An = 4.220. 
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Table 1 



The observed and calculated vibrational levels are presented in Table 1 where our 
results are also compared with those calculated in the boson-reahzation model 
The AEf-ai arc the differences between the observed and the calculated values. The 
standard deviation (SD) of this fit is 6.724 cm~^, while that of the boson-realization 
model and the corresponding gr-deformed model is 8.198 cm~^ and 7.157 cm~^ 1^°], 
respectively. The present model applied to H2O gives a slight improvement in the fit 
to its vibrational energy levels. 

2.2 Application to ASH3 

To our knowledge, vibrations of a molecule XY3 have not been studied via algebraic 
methods yet. As an example, we now apply Eq. (1.5) to the vibrational spectrum 
of arsine (AsHs). The calculation of energy levels will be greatly simplified if the 
symmetrized bases are used. For the considered system, these bases can be obtained 
by the method of Ref . [7] . This method is particularly useful when describing vibrations 
of polyatomic molecules and for high overtones. 

Recently observed data (Eobs) of AsHa with their normal labels {1^11^21^3^ are 
taken from the compilation of Refs. [9,15]. Since there are 18 experimental data for 
this molecule it is not possible to determine all coupling constants in Eq. (1.5) by 
a least-squares optimization. It is natural to choose the following parameters in Eq. 
(1.5) nonvanishing: Ug, Xg, cj^, x^, rji, ri2, rj^, rj^, rj^, and rjQ. In other words, we neglect 
the terms with fjt {7 <t < 11). We have found out that other choice does not provide 
improvement in our fits. For comparison, we have made the first fit (Fit a) where the 
model is in the limit of A^^ — > 00 and A^b — > 00. The obtained parameters and the 
calculated differences (AE^^i) with the corresponding observed data (Eobs) in this fit 
are listed in Table 2 and 3. In another fit (Fit b) we take the boson numbers Ns and 
Nf, to be 62 and 42, respectively. The results in Fit b are also listed in the related 
Tables. 

Table 2 
Table 3 



6 



The SD in Fit a and b is 10.201 cm~^ and 5.827 cm~\ respectively. The calculated 
results for this molecule in Fit b can be compared with those by Lukka et al '^1 . However, 
their vibrational Hamiltonian was expressed in terms of curvilinear internal coordinates, 
and the Morse oscillator basis functions and the results of ah initio calculation were 
used there. 

3. Conclusion and discussion 

We have proposed Fermi resonance- algebraic models for the vibrations of a bent 
molecule XY2 and a molecule XY3 by U(2) algebras for the description of both stretch- 
ing and bending modes, where Fermi resonances between the stretch and the bend are 
considered. In the limit, the models reduce to the recently presented boson-realization 
models. As an example, they have been applied to vibrational spectra of the water 
molecule and arsine (AsHa). To some extend, they have provided improvement in the 
standard deviations in the corresponding boson-realization models for the same ob- 
served energy levels. It needs more sample calculations to judge which approach is 
better. 

The present calculations and others [^^'^'I'^'i^l demonstrate that an algebraic ap- 
proach can be used as an effective method for describing molecular vibrations with 
good precision. In this approach the eigenvalues and the related wave functions are 
obtained through matrix diagonalization. Hence, the required computing time is short, 
and some other physical properties, such as transition intensities I^'^^l, are allowed to 
calculate. It provides the possibility of using an algebraic model as a numerically effi- 
cient, phenomenological theory for studying molecular spectra, especially when no ah 
initio calculation is available. 
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TABLE 1. Observed and calculated vibrational energy levels of H2O ( cm ^) 





^ohs 




A TP 


(l/lZ^2'^3) 


^obs 


A pIlO] 

cat 


A TP 


(010) 


1594.7498 


-2.668 


-4.733 


(001) 


3755.93 


0.170 


-0.522 


(020) 


3151.63 


0.016 


-4.136 


(Oil) 


5331.269 


-5.434 


-5.562 


(100) 


3657.053 


3.549 


1.997 


(021) 


6871.51 


-3.398 


-4.325 


(030) 


4666.793 


5.464 


0.110 


(101) 


7249.81 


2.611 


0.409 


(110) 


5234.977 


1.684 


0.694 


(031) 


8373.853 


4.051 


2.443 


(040) 


6134.03 


9.088 


4.828 


(111) 


8807 


-1.858 


-0.673 


(120) 


6775.1 


1.817 


-0.591 


(041) 


9833.584 


13.643 


12.564 


(200) 


7201.54 


1.302 


0.483 


(121) 


10328.731 


-0.616 


-0.006 


(002) 


7445.07 


3.063 


1.453 


(201) 


10613.355 


0.766 


-0.809 


(050) 


7552 


11.254 


10.415 


(003) 


11032.406 


5.575 


3.577 


(130) 


8273.976 


0.594 


-1.912 


(131) 


11813.19 


2.636 


1.919 


(210) 


8761.582 


0.248 


1.828 


(211) 


12151.26 


-3.666 


0.371 


(060) 


8890 


-17.089 


-12.506 


(013) 


12565 


-9.461 


-6.151 


(012) 


9000.136 


-5.600 


-4.650 


(141) 


13256 


4.292 


4.464 


(220) 


10284.367 


-0.405 


-0.714 


(221) 


13652.656 


-3.450 


-0.116 


(022) 


10524.3 


-0.950 


-0.991 


(301) 


13830.938 


0.384 


-1.448 


(300) 


10599.686 


-3.117 


-3.002 


(023) 


14066.194 


-9.806 


-7.824 


(102) 


10868.876 


5.583 


1.810 


(103) 


14318.813 


6.331 


1.985 


(310) 


12139.2 


-6.313 


-1.508 


(151) 


14640 


-10.977 


-5.968 


(112) 


12407.64 


3.124 


3.469 


(231) 


15119.029 


-3.536 


-2.852 


(042) 


13448 


7.777 


4.483 


(311) 


15347.956 


-4.276 


1.217 


(320) 


13642.202 


-7.521 


-4.024 


(113) 


15832.765 


-7.239 


-1.381 


(202) 


13828.277 


2.052 


0.344 


(321) 


16821.635 


-2.039 


2.715 


(122) 


13910.896 


6.188 


-0.130 


(203) 


16898.842 


3.234 


-1.570 


(400) 


14221.161 


-1.522 


-1.806 


(123) 


17312.539 


1.834 


-1.041 


(004) 


14536.87 


9.844 


6.393 


(401) 


17495.528 


0.756 


-0.563 


(330) 


15107 


-14.579 


-12.251 


(331) 


18265.82 


-3.654 


-2.516 


(212) 


15344.503 


12.518 


10.038 


(213) 


18393.314 


-1.400 


0.480 


(410) 


15742.795 


0.179 


4.839 


(411) 


18989.961 


-12.986 


-0.621 


(222) 


16825.23 


-16.629 


-18.041 


(421) 


19720 


-0.801 


-0.492 


(302) 


16898.4 


1.469 


-3.711 


(303) 


19781.105 


14.596 


8.509 


(420) 


17227.7 


11.514 


0.715 


(501) 


20543.137 


-7.555 


-2.973 


(500) 


17458.354 


0.824 


0.139 


(313) 


21221.828 


-6.455 


-5.430 


(104) 


17748.073 


-0.279 


-0.504 


(232) 


18320 


25.584 


29.589 


(312) 


18392.974 


-1.070 


0.979 


(412) 


21221.569 


-5.258 


-5.451 












SD 


8.198 


6.724 
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TABLE 2. Parameters (cm ^) for AsHa obtained by the least square fitting 



CUs Xs CUb Xb Tji 

Fit a 2039.921 -80.580 959.021 64.320 -24.775 
Fit b 2066.972 -28.376 962.647 68.865 -18.639 



m m m % ve 

Fit a -21.932 -53.616 -5.996 -33.712 1.823 
Fit b -27.283 -55.482 -5.850 -29.862 11.970 



TABLE 3. Observed and calculated vibrational energy levels of AsHs (cm ^) 



(l^ll/2t'3V4-') 


f[9,15] 
^obs 








p[9,15] 
^obs 






(OIOOqO) 


906.752 


-8.425 


-1.328 


(0000l±i) 


999.225 


8.253 


9.295 


(020000) 


1806.149 


3.950 


1.853 


(0100l±i) 


1904.115 


-9.368 


-0.427 


(000020) 


1990.998 


-2.224 


1.602 


(00002±2) 


2003.483 


-10.736 


-9.707 


(IOOOqO) 


2115.164 


-0.720 


0.932 


(001^0°) 


2126.423 


2.698 


-1.537 


(IIOOqO) 


3013 


1.704 


3.953 


(00lH±i) 


3102 


-0.615 


3.396 


(20000°) 


4166.772 


-0.182 


3.033 


(lOl^QO) 


4167.935 


-8.389 


-2.245 


(SOO^OO) 


6136.316 


2.401 


-4.667 


(20110°) 


6136.310 


9.419 


1.951 


(OO3O0O) 


6276 


-3.991 


-2.321 


(10220°) 


6295 


2.410 


2.238 


(40000°) 


8028.977 


-2.275 


-1.389 


(301^0°) 


8028.969 


1.398 


0.713 



SD 10.201 5.827 
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